Method of estimating a signal-to-interferenceratio (SINR)

ABSTRACT

In the method of estimating a signal-to-interference+noise ratio (SINR), an initial SINR estimate is generated based on a mean of a plurality of samples and a sample variance estimate of the plurality of samples. Then, the initial SINR estimate is scaled and translated. Furthermore, SINR estimates based on pilot symbols and data symbols, respectively, are combined to form a composite estimate.

BACKGROUND OF THE INVENTION

[0001] 1. Field of the Invention

[0002] The present invention relates to wireless communication, and more particularly, a method of estimating a signal-to-interference+noise ratio.

[0003] 2. Description of Related Art

[0004] Signal-to-Interference+Noise Ratio (SINR) is an important metric of communication link quality. SINR estimation is of particular importance for wireless data systems where resources are shared dynamically amongst users. Some applications of SINR estimates are: a) Power Control in CDMA Systems: the receiver estimates the SINR, compares it to a target and commands the transmitter to increase/decrease its transmitted power; and b) Rate Adaptation: the information bit-rate assigned to a user can be dynamically varied based on its link quality and the system load. While such adaptation has limited use in voice systems, it is extremely useful for wireless data systems. Consequently, inaccurate SINR estimates can severely degrade performance and resource utilization.

SUMMARY OF THE INVENTION

[0005] In the method of estimating a signal-to-interference+noise ratio (SINR) according to the present invention, a conventionally known SINR estimate is scaled and translated to produce a SINR estimate that has a substantially reduced mean square error. The SINR estimator according to the present invention operates on either pilot symbol samples or data symbol samples. Further according to the present invention, a pilot symbol based SINR estimate and a data symbol based SINR estimate are weighted and combined to generate a composite SINR estimate having a reduced mean square error as compared to either of the individual estimators.

BRIEF DESCRIPTION OF THE DRAWINGS

[0006] The present invention will become more fully understood from the detailed description given herein below and the accompanying drawings, which are given by way of illustration only, and thus are not limitative of the present invention, and wherein:

[0007]FIG. 1 illustrates pilot and data symbols multiplexed in a time slot;

[0008]FIG. 2 illustrates the gain from the SINR estimator according to an embodiment of the present invention as compared to a conventional estimate based on pilot symbols only for different smoothing factors.

[0009]FIG. 3 illustrates optimal weights assigned to a pilot-symbol based estimator and a data-symbol based estimator in a combined estimators; and

[0010]FIG. 4 illustrates the performance of the combined estimator using the estimated weights.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0011] In describing the method of estimating the signal-to-inteference+noise ratio according to the present invention, only Binary Phase Shift Keying (BPSK) modulation is considered although the methods and related analysis can be extended to other signaling schemes. Noise and interference are modeled together as additive white Gaussian noise (AWGN), but as will be appreciated from the following disclosure, this should not limit the application of the method according to the present invention. Transmission is organized into fixed duration timeslots, each containing pilot and data symbols as shown in FIG. 1. The channel attenuation and phase shift, assumed fixed over a timeslot, are treated as unknown constants rather than as random variables (no fading assumption).

[0012] First, to better understand the method according to the present invention, the conventional SINR estimation will be described. Typically, the received signal corresponding to the jth transmitted symbol (pilot or data) in the kth time slot is defined as

Y _(kj) =a _(kj)μ_(k)+ε_(kj) j=1, 2, . . . , N,  (1)

[0013] where μ_(k) represents the received signal amplitude (product of transmitted amplitude and channel gain), ε_(kj) is a random variable that represents the noise+interference, a_(kj) represents the symbol-value, and N is the number of samples (pilot or data). Information symbols could be +1 or −1 (in BPSK), while it is assumed (without any loss of generality) that pilot symbols are always +1. It is also assumed that the distribution that characterizes the noise+interference is Gaussian with zero mean and variance σ² The SINR in the kth time slot is defined as: $\begin{matrix} {\Theta_{k} = \frac{\mu_{k}^{2}}{\sigma^{2}}} & (2) \end{matrix}$

[0014] and is the parameter to be estimated.

[0015] The groups of N sample points (data or pilot) could correspond to a time slot in CDMA systems or frames in TDMA systems. A well-known pilot-symbol sample based estimator of SINR is computed as the ratio of the square of the sample mean of the received to the pilot-symbol sample Y (based on N sample points in a group) to the corresponding sample variance is a reasonable estimator of the SINR. Estimators based on this ratio are called Squared Mean By Variance or SMV estimators. Different SMV estimators have been studied in the literature and they only differ in the normalization constant used to compute the sample variance.

[0016] For the case where the {Y_(kj)} values correspond to pilot symbols, define the sample mean and unbiased sample variance for the kth time slot as $\begin{matrix} {{\overset{\_}{Y}}_{k} = {\frac{1}{N}{\sum\limits_{j = 1}^{N}Y_{kj}}}} & (3) \\ {S_{k}^{2} = {\frac{1}{N - 1}{\sum\limits_{j = 1}^{N}{\left( {Y_{kj} - {\overset{\_}{Y}}_{k}} \right)^{2}.}}}} & (4) \end{matrix}$

[0017] Then, ${\hat{\Theta}}_{k} = {\left( {\overset{\_}{Y}}_{k} \right)^{2}/S_{k}^{2}}$

[0018] is one possible SMV estimator (commonly used). The Maximum Likelihood (ML) estimator of the SINR is also an SMV estimator where a biased sample variance estimate is used in the ratio (normalization is by N rather than N-1). The signal power, μ_(k) ², varies due to channel fading. However, the noise variance, σ², changes very slowly with time, typically with the addition (departure) of a call; therefore, one can improve the overall quality of the SINR estimate by using a longer-term moving average estimate of the sample variance. This makes the “effective” number of samples used in the sample variance estimate larger and therefore more accurate. One simple method to accomplish this is through exponential smoothing of a set of sample variance estimates (called Exponentially Weighted Moving Average or EWMA). The smoothed sample variance estimate through the kth time slot would be given by: $\begin{matrix} \begin{matrix} {{\hat{\sigma}}_{k}^{2} = {{\left( {1 - r} \right){\hat{\sigma}}_{k - 1}^{2}} + {r\quad S_{k}^{2}}}} & {{k \geq 1},} \end{matrix} & (5) \end{matrix}$

[0019] where r is the smoothing factor determined according to desired design parameters and 0<r≦1. The SINR estimate at the end of k time slots then becomes: $\begin{matrix} {{\hat{\Theta}}_{k} = {\frac{\left( {\overset{\_}{Y}}_{k} \right)^{2}}{{\hat{\sigma}}_{k}^{2}}.}} & (6) \end{matrix}$

[0020] The accuracy (mean and variance) of SMV estimators based on using a single group of pilots with N sample points is known in the art. However, aside from simulation results, the accuracy of SMV estimators that utilize EWMA for sample variance estimation is not known.

[0021] The inventors have analytically analyzed the accuracy of SMV estimators that utilize EWMA for the sample variance estimate and have developed accurate approximations for the bias, variance and MSE of such estimators. From this analysis, the inventors developed a characterization of the bias and variance, and hence the MSE, of such estimators. The inventors, based on their analysis, developed a scaling and translation methodology for the SINR estimators that substantially reduces the MSE. The scaled and translated EWMA SINR estimator according to the present invention is: $\begin{matrix} {{{\overset{\sim}{\Theta}}_{k} = {{\frac{\eta - 2}{\eta}{\hat{\Theta}}_{k}} - \frac{1}{N}}},} & (7) \end{matrix}$

[0022] where $\left( {- \frac{1}{N}} \right)$

[0023] is the translation and η is the scaling factor defined as: $\begin{matrix} {\eta = {\frac{\left( {N - 1} \right)\left( {2 - r} \right)}{r}.}} & (8) \end{matrix}$

[0024] Note that this scaling factor η depends only on parameters that are fixed and known at the receiver and do not require any separate estimation procedure.

[0025] The gains due to scaling and translation (measured as the reduction in mean-squared error) in using {tilde over (Θ)}_(k) instead of {circumflex over (Θ)}_(k) can be substantial especially for larger values of the smoothing factor r as seen in FIG. 2. FIG. 2 illustrates the gain (measured as the reduction in MSE) based on pilot symbols only for different smoothing factors r, and show that even for values of r around 0.5-0.6, the gain is around 30-40%.

[0026] SMV estimators that use data symbols are also known in the art. A commonly used SMV estimator based on data symbols, called the non-coherent estimator, is one where the sequence {Y_(kj)} is replaced by the sequence of its absolute values i.e. {Z_(kj)=|Y_(kj)|}. The SINR estimate for the kth time slot is the ratio of the sample mean and sample variance of the sequence {Z_(kj)}. Smoothing of the sample variance of the {Z_(kj)} sequence via the EWMA approach may also be used to improve accuracy. While the scaling and translation approach developed in this invention can be extended to any SMV estimator based on data symbols, including the non-coherent estimator, the preferred approach is outlined in concurrently filed Application No. ______ entitled “Method of Estimating a Signal-to-Interference+Noise Estimate Using Data Samples” by the inventors of the subject application.

[0027] The accuracy of the SINR estimate can be further enhanced by simultaneously utilizing estimates from pilot symbols as well as data symbols in the estimation process. As mentioned before pilot symbols have the advantage of being known at the receiver but are relatively fewer in number. Data symbols are more plentiful than pilot symbols but they are unknown at the receiver.

[0028] First, a method for the general problem of optimally linearly combining two estimates of an unknown parameter will be described. With that foundation, the combining method will be applied to SINR estimation with pilot and data symbols.

[0029] Suppose there are two noisy estimates, {circumflex over (Θ)}₁ and {circumflex over (Θ)}₂, respectively, of a parameter Θ. Furthermore let the biases of the two estimates be B₁=E[{circumflex over (Θ)}₁]−Θ and B₂=E[{circumflex over (Θ)}₂]−Θ, respectively, and the variances be τ₁ ² and τ₂ ², respectively. The MSE of {circumflex over (Θ)}₁ is MSE₁=B₁ ²+τ₁ ², i=1,2. Consider the class of estimators formed by an arbitrary linear combination of {circumflex over (Θ)}₁ and {circumflex over (Θ)}₂

{circumflex over (Θ)}=a ₁{circumflex over (Θ)}₁ +a ₂{circumflex over (Θ)}₂.  (9)

[0030] Minimum mean squared error weights a₁ and a₂ are obtained by minimizing E({circumflex over (Θ)}−{circumflex over (Θ)})². The following are the minimum mean squared error weights $\begin{matrix} {\alpha_{1}^{o\quad p\quad t} = {\frac{\left( {\Theta + B_{1}} \right)\tau_{2}^{2}}{{\tau_{1}^{2}\left( {\Theta + B_{2}} \right)}^{2} + {\tau_{2}^{2}\left( {\Theta + B_{1}} \right)}^{2} + {\tau_{1}^{2}\tau_{2}^{2}}}\Theta}} & (10) \\ {\alpha_{2}^{o\quad p\quad t} = {\frac{\left( {\Theta + B_{2}} \right)\tau_{1}^{2}}{{\tau_{1}^{2}\left( {\Theta + B_{2}} \right)}^{2} + {\tau_{2}^{2}\left( {\Theta + B_{1}} \right)}^{2} + {\tau_{1}^{2}\tau_{2}^{2}}}{\Theta.}}} & (11) \end{matrix}$

[0031] The corresponding minimum MSE is $\begin{matrix} {{M\quad S\quad E^{o\quad p\quad t}} = {\frac{\Theta^{2}}{1 + {\left( {\Theta + B_{1}} \right)^{2}/\tau_{1}^{2}} + {\left( {\Theta + B_{2}} \right)^{2}/\tau_{2}^{2}}}.}} & (12) \end{matrix}$

[0032] Unfortunately, the optimum weights to use depend on the value of Θ. Later in this section, we shall describe how the weights can be estimated in practice.

[0033] When two unbiased estimators are being combined (i.e. B₁=0 and B₂=0) and it is desired that the combined estimators also be unbiased, then the weights need to be constrained to sum to unity. In this case the optimal weight on {circumflex over (Θ)}₁, can be shown to be $\begin{matrix} {\alpha = \frac{{M\quad S\quad E_{2}} - {B_{1}B_{2}}}{{M\quad S\quad E_{1}} + {M\quad S\quad E_{2}} - {2B_{1}B_{2}}}} & (13) \end{matrix}$

[0034] and the optimal weight on {circumflex over (Θ)}₂ is thus 1−a=(MSE₁−B₁B₂)/(MSE₁+MSE₂−2B₁B₂). The MSE of the optimal constrained estimator is $\begin{matrix} {{M\quad S\quad E^{o\quad p\quad t}} = {\frac{{M\quad S\quad E_{1}M\quad S\quad E_{2}} - {B_{1}^{2}B_{2}}}{{M\quad S\quad E_{1}} + {M\quad S\quad E_{2}} - {2B_{1}B_{2}}}.}} & (14) \end{matrix}$

[0035] When both estimators are unbiased Equation (13) reduces to the known result for optimally combining (linearly) two independent unbiased estimates.

[0036] The preceding results will now be applied to the problem of SINR estimation. Consider two estimators of the SINR, Θ_(k). For the purpose of illustration we select the smoothed, scaled and translated pilot-sample based estimator, say {tilde over (Θ)}_(k) ^(P), and the smoothed, scaled and translated data-sample based estimator (e.g. the non-coherent estimator), say {tilde over (Θ)}_(k) ^(D), as the two estimators to be combined. The inventors have shown that the biases and variances of these estimators are approximated as follows $\begin{matrix} {{B_{k}^{P} = 0}{\left( \tau_{k}^{P} \right)^{2} = {{\frac{2\left( {\eta_{P} - 2} \right)}{N_{P}^{2}}\left\lbrack {\frac{\left( {1 + {N_{P}\Theta_{k}}} \right)^{2}}{\left( {\eta_{P} - 2} \right)\left( {\eta_{P} - 4} \right)} + \frac{1 + {2N_{P}\Theta_{k}}}{\eta_{P} - 4}} \right\rbrack}\left( {\eta_{P} > 4} \right)}}{B_{k}^{D} = {{\Theta_{z,k} - {\Theta_{k}\left( \tau_{k}^{D} \right)}^{2}} = {{\frac{2\left( {\eta_{D} - 2} \right)}{N_{D}^{2}}\left\lbrack {\frac{\left( {1 + {N_{D}\Theta_{Z,i}}} \right)^{2}}{\left( {\eta_{D} - 2} \right)\left( {\eta_{D} - 4} \right)} + \frac{1 + {2N_{D}\Theta_{Z,i}}}{\left( {\eta_{D} - 4} \right)}} \right\rbrack}{\left( {\eta_{D} > 4} \right).}}}}} & (15) \end{matrix}$

[0037] In Equation (15), the superscript and subscript P denotes the pilot sample based variables, the superscript and subscript D denotes the data sample based variables, N_(D) and N_(P) are the number of data and pilot symbols respectively per time slot, η_(P)=(N_(P)−1)(2−r)/r, η_(D)=(N_(D)−1)(2−r)/r (where r is the EWMA smoothing factor) and the quantity where ${\Theta_{z,k} = \frac{\mu_{Z\quad k}^{2}}{\sigma_{Z,k}^{2}}},{w\quad h\quad e\quad r\quad e}$ $\mu_{Z\quad k} = {{\sqrt{\frac{2}{\pi}}{\exp \left( {- \frac{\Theta_{k}}{2}} \right)}} + {\frac{\sqrt{\Theta_{k}}}{2}\left\lbrack {{{erfc}\left( {- \sqrt{\frac{\Theta_{k}}{2}}} \right)} - {{erfc}\left( \sqrt{\frac{\Theta_{k}}{2}} \right)}} \right\rbrack}}$ σ_(Z  k)² = Θ_(k) + 1 − μ_(Z  k)²

[0038] (erfc in the above equation is the well-know complementary error function). Using Equation (15) with either Equation (10) and Equation (11) or Equation (13) gives the optimal combined estimators of {tilde over (Θ)}_(k) ^(P) and {tilde over (Θ)}_(k) ^(D). Since the biases and variances, and therefore the weights themselves depend on Θ_(k), which is unknown, one can only estimate the weights in practice. An estimate of the weights can be obtained by replacing Θ_(k) by {tilde over (Θ)}_(k) ^(P) in the corresponding expressions. One could, alternatively, replace Θ_(k) by {tilde over (Θ)}_(k) ^(D) to obtain estimates of these quantities. However, since {tilde over (Θ)}_(k) ^(P) is an unbiased estimate of Θ_(k), the accuracy of the composite estimate is better, when {tilde over (Θ)}_(k) ^(P) is used. Note that the above technique for combining may be used with any two estimators as long as the appropriate bias and variance expressions are used in the weight computation. For example, one can combine unsmoothed, scaled and translated pilot and data-sample based estimators using the same approach.

[0039] Let C1 denote the combined estimator of {tilde over (Θ)}_(k) ^(P) and {tilde over (Θ)}_(k) ^(D) that has constrained weights and C2 denote the combined estimator with unconstrained weights. The optimal weights used when {tilde over (Θ)}_(k) ^(P) and {tilde over (Θ)}_(k) ^(D) are combined for C1 and C2 are shown in FIG. 3, as a function of Θ_(k). In the figure-legend BCSV (BCSV-Z) stands for the smoothed, scaled and translated pilot (data) sample based estimator. The root mean squared error (RMSE) values for C1 and C2 when estimated weights are used are shown in FIG. 4. For comparison, the corresponding RMSE values for each of the component estimators, BCSV and BCSV-Z, are also shown. (The label “A” in the figure merely means that numerical results were obtained from analytic formulas rather than simulations). Despite the error in estimating the weights, both combined estimators outperform the individual estimators for the entire range of SINR considered.

[0040] The invention being thus described, it will be obvious that the same may be varied in many ways. Such variations are not to be regarded as a departure from the spirit and scope of the invention, and all such modifications as would be obvious to one skilled in the art are intended to be included within the scope of the following claims. 

We claim:
 1. A method of estimating a signal-to-inteference+noise ratio (SINR), comprising: generating an initial SINR estimate based on a mean of a plurality of samples and a sample variance estimate of the plurality of samples; scaling the initial SINR estimate; and translating the scaled SINR estimate.
 2. The method of claim 1, wherein the generating step generates the initial SINR estimate based on at least two sample variance estimates.
 3. The method of claim 2, wherein the generating step generates a smoothed sampled variance estimate based on the sample variance estimates, and generates the initial estimate based on the smoothed sample variance estimate.
 4. The method of claim 3, wherein the generating step generates the smoothed sampled variance estimate by combining a current sample variance estimate and a previous sample variance estimate based on a smoothing factor.
 5. The method of claim 4, wherein the scaling step scales the initial SINR estimate based on the smoothing factor.
 6. The method of claim 5, wherein the generating step generates the smoothed sample variance estimate according to the following expression: σ̂_(k)² = (1 − r)σ̂_(k − 1)² + r  S_(k)²

where {circumflex over (σ)}_(k) ²=current smoothed sample variance estimate, {circumflex over (σ)}_(k−1) ²=previous smoothed sample variance estimate, S_(k) ²=current sample variance estimate, and r=smoothing factor.
 7. The method of claim 5, wherein the scaling step scales the initial SINR estimate based on the following expression: $\overset{\sim}{\Theta} = {\frac{\eta - 2}{\eta}\hat{\Theta}}$

where {tilde over (Θ)}=scaled SINR estimate, {circumflex over (Θ)}=initial SINR estimate, and $\eta = {{\frac{\left( {N - 1} \right)\left( {2 - r} \right)}{r}\quad {where}\quad N} = {a\quad {number}\quad {of}\quad {the}\quad {plurality}\quad {of}\quad {{samples}.}}}$

where N=a number of the plurality of samples.
 8. The method of claim 7, wherein the translating step translates the scaled SINR estimate based on the number of the plurality of samples.
 9. The method of claim 8, wherein the translating step translates the scaled SINR estimate by adding (−1/N), where N is the number of the plurality of samples.
 10. The method of claim 4, wherein the scaling step scales the initial SINR estimate based on a number N of the plurality of samples.
 11. The method of claim 1, wherein the translating step translates the scaled SINR estimate based on the number of the plurality of samples.
 12. The method of claim 11, wherein the translating step translates the scaled SINR estimate by adding (−1/N), where N is the number of the plurality of samples.
 13. The method of claim 1, wherein the plurality of samples are pilot symbol samples.
 14. The method of claim 1, wherein the plurality of samples are data symbol samples.
 15. A method of estimating a signal-to-inteference+noise ratio (SINR), comprising: generating a first initial SINR estimate based on a mean of a plurality of pilot symbol samples and an sample variance estimate of the plurality of pilot symbol samples; scaling the first initial SINR estimate; translating the first scaled SINR estimate; generating a second initial SINR estimate based on a mean of a plurality of data symbol samples and an sample variance estimate of the plurality of data symbol samples; scaling the second initial SINR estimate; translating the second scaled SINR estimate; and combining the first and second scaled estimates to produce a composite SINR estimate.
 16. The method of claim 15, wherein the combining step weights the first and second scaled estimates and combines the weighted first and second scaled estimates to produce the composite SINR estimate.
 17. The method of claim 15, wherein the weights are determined based on a bias in the first and second scaled estimates, respectively.
 18. A method of estimating a signal-to-inteference+noise ratio (SINR), comprising: generating a first SINR estimate based on received pilot symbol samples; generating a second SINR estimate based on received data symbol samples; and combining the first and second SINR estimates to produce a composite SINR estimate.
 19. The method of claim 18, wherein the combining step comprises: weighting the first SINR estimate according to a first weight; weighting the second SINR estimate according to a second weight; and combining the first and second weighted SINR estimates.
 20. The method of claim 19, wherein the first and second weights sum to unity.
 21. The method of claim 19, wherein the first and second weights do not sum to unity.
 22. The method of claim 19, wherein the first and second weights depend on an initial estimate of the SINR.
 23. The method of claim 19, wherein the first weight depends on a bias and variance of the first SINR estimate and the second weight depends on a bias and variance of the second SINR estimate.
 24. The method of claim 19, wherein the generating the first SINR estimate step includes, generating a first initial SINR estimate based on a mean of a plurality of pilot symbol samples and a sample variance estimate of the plurality of pilot symbol samples, scaling the first initial SINR estimate, and translating the first scaled SINR estimate; and the generating the second SINR estimate step includes, generating a second initial SINR estimate based on a mean of a plurality of data symbol samples and a sample variance estimate of the plurality of data symbol samples, scaling the second initial SINR estimate, and translating the second scaled SINR estimate. 